One particle in a box: the simplest model for a Fermi gas in the unitary limit 
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We consider a single quantum particle in a spherical box interacting with a fixed scatterer at the 
center, to construct a model of a degenerate atomic Fermi gas close to a Feshbach resonance. One 
of the key predictions of the model is the existence of two branches for the macroscopic state of the 
gas, as a function of the magnetic field controlling the value of the scattering length. This model 
is able to draw a qualitative picture of all the different features recently observed in a degenerate 
atomic Fermi gas close to the resonance, even in the unitary limit. 
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Recently several experiments have produced two spin- 
component degenerate Fermi gases in the unitary limit, 
that is in a limit where the scattering length a of two 
atoms with different spin components is much larger than 
the mean interparticle separation P, 0, J,, .4, 5j. Simi- 
lar experiments have also been performed with bosonic 
atoms 0, 01 Hi 13 ■ This was made possible by a tuning 
of the scattering length virtually from —00 to -l-oo us- 
ing a Feshbach resonance driven by an external uniform 
ma gnet ic field, as first demonstrated on bosonic atoms 

This regime constitutes a theoretical challenge. Some 
theories, based on the Hartree-Fock mean field approx- 
imation for the normal gas or on the BCS mean field 
calculation for the superfluid phase, rely on the small 
parameter kpa^ where the Fermi momentum of the gas 
is conventionally related to the mean total density p by 
the non-interacting case formula: 
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Such approaches are not quantitative in the limit kF\a\ 
+QO. More sophisticated approaches have been devel- 
oped to extend the accuracy of the theory to the unitary 
regime Elll El El El . The difliculty comes from 
the fact that there exists no obvious small parameter for 
the theory in the unitary limit, at least in the degenerate 
regime of a temperature T much smaller than the Fermi 
temperature Tp , which is the regime considered here and 
also the present experimental situation. Note that, in the 
classical regime T ^ Tp, one recovers a small parameter 
/o^/'^l/l <^C 1 since the typical scattering amphtude / for 
\a\ — > -l-oo, is on the order of A, where A is the thermal 
de Broglie wavelength 0, . 

The scope of the present work is to give the simplest 
possible physical picture of a two spin-component Fermi 
gas for arbitrary values of the scattering length a, essen- 
tially at zero temperature. The model is not made to give 
quantitative predictions, but, as we shall see, it qualita- 
tively reproduces the experimental observations, and it 
will provide, we hope, useful guidelines for experiments 
to come. 



The model: We consider a spatially homogeneous gas 
of N/2 fermions of spin -1-1/2 and N/2 fermions of spin 
— 1/2, each particle having a mass m. The interaction 
of a given spin -1-1/2 particle with the N/2 spin —1/2 
particles is modelized (i) by the interaction of a fictitious 
particle of mass equal to the reduced mass m/2, with a 
fixed scatterer at the center of a spherical box of radius R, 
and (ii) with the boundary conditions that the wavefunc- 
tion 4){r) of the fictitious particle vanishes on the surface 
of the box. (i) : The interaction of the fictitious particle 
with the scatterer represents the interaction of the given 
spin -1-1/2 particle with its nearest spin —1/2 neighbour: 
in the model, these two opposite spin particles are con- 
sidered in their center of mass frame and in the singlet 
spin state, (ii): The boundary condition mimics the in- 
teraction effect of the N /2 — 1 other spin —1/2 particles 
and the Fermi statistical effect of the remaining N/2 — 1 
spin -1-1/2 particles, see Figure ^ A similar model was 
very recently put forward for bosons, with the difference 
that the box is replaced by a harmonic potential [T9|. 

In the absence of interactions between the fermions, 
the total energy of the gas is given by the known ideal 
Fermi gas formula in the thermodynamic limit: 
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where the Fermi energy ep = fi'^kp/2m is related to the 
mean density through The gas energy E is also re- 
lated to the energy e of the fictitious particle in the spher- 
ical box by 



E = -Ne 
2 



(3) 



For the ideal Fermi gas, there is no scatterer in the box 
so that the ground state value of e is 
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This establishes the link between the radius R and the 
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FIG. 1: Model used: a fictitious particle corresponding to the 
relative motion of a fermion with spin +1/2 and the nearest 
fermion with spin —1/2 is scattered by the fixed delta Fermi 
pseudo-potential [2^ in the center of a spherical box of radius 
R with absorbing boundary conditions. The box mimics the 
interaction effect of the N/2 — 1 remaining spin —1/2 atoms 
and the Fermi statistical effect of the remaining N/2 — 1 spin 
+1/2 fermions. The coupling constant g is 4iYh^a/m, where 
a is the scattering length of two fermions with opposite spin, 
g5ieg{r) = gS{f)dr{r.) is the Fermi pseudo-potential which 
imposes the contact condition ||5J on the two-body wave func- 
tion 



mean density of the gas [23 : 
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In presence of interactions, the scatterer in the model 
is the Fermi delta pseudo-potential with a coupling con- 
stant g = ATrh^a/m where a is the s-wave scattering 
leng th of two fermions with opposite spin components 
[21I I22 . |23j I2JI ■ The choice of such a zero-range poten- 
tial is allowed in the regime 



pri < 1 
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where r,, is the effective range of the true interaction po- 
tential, and which allows to consider these systems as 
gases rather than liquids even in the unitary limit. The 
two-body scattering amplitude derived from the pseudo- 
potential. 
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is indeed a good approximation of the scattering ampli- 
tude in simple models of a Feshbach resonance: 
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when fej^lfel <C 1 since k is at most on the order oi kp 
at T <C Tp. Te was calculated in and one finds that 
condition © is very well satisfied in present experiments 
on Li6 ^151. 

We will use the fact that the pseudo-potential is equiv- 
alent to re placing the interaction potential by the contact 
condition '2lll2al2l 
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where r is the distance to the origin. Out of the ori- 
gin in the box, the wavefunction then solves the free 
Schrodinger equation: 
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Restricting to the s-wave, which is the only partial wave 
in which the pseudo-potential scatters, is rotationally 
symmetric and we set 0(r) — u{r)/r. For positive ener- 
gies e = Ti^fc^/m, where fc > 0, one then finds 



u{r) oc sin[A;(r — K)] , 
and the contact condition Eq.lO imposes: 

tan kR =^ ka 
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For negative energies e — —h^K^/m, where k > 0, one 
finds 



u{r) oc sinh[fc(r — R)] 
where Eq.(|5J) imposes: 

tanh kR — na . 



(13) 
(14) 




FIG. 2: Total energy of the gas in units of the Fermi energy, 
as function of —l/kpa. Only the first two branches of the 
model are represented. 

Discussion: The energy spectrum of the previous model 
allows to calculate the energy of the gas according to 



Eq.Q. The result is shown in Figure |2] as function of 
— 1/kpa, restricting to the two lowest energy branches. 
The choice of the parameter —1/kpa is inspired by the 
way the experimental results on ^Li are usually pre- 
sented: in this way, the right part of the Figure |2] corre- 
sponds to magnetic fields above the Feshbach resonance, 
and the left part to magnetic fields below the resonance. 
The existence of several energy branches for the macro- 
scopic state of the gas is a first crucial prediction of the 
model. The ground branch connects two clearly identi- 
fied regimes: 

• the first regime kpa — s- 0^ (on the right hand side 
of Figure 2) is a weakly attractive Fermi gas which 
corresponds at zero temperature to the BCS phase 
in a more complete treatment. 

• the second regime kpa ^ (on the left hand side 
of Figure 2) corresponds to a dilute gas of dimers; 
each dimer corresponds to the bound state of the 
two-body problem in free space since the size of 
the box R, on the order of the mean interatomic 
distance, is here much larger than the spatial ex- 
tension '--^ a of the dimer. The energy of a dimer 
is —h^ /ma? which explains the quadratic drop of 
the ground energy on the left part of the reso- 
nance. In a more complete treatment, one would 
find at r = that these dimers form a Bose- 
Einstein condensate since they are bosons in the 
limit pi/^a < 1. 

The upper branch in Figure |21 is also easy to identify in 
the weakly interacting limit kpa — > 0^, where it corre- 
sponds to a weakly repulsive Fermi gas. It clearly cor- 
responds to a metastable state of the gas: a three body 
collision between atoms will form a dimer plus an ex- 
tra atom carrying away the binding energy, a mechanism 
that will depopulate the upper branch and populate the 
ground branch. This constitutes a first experimental way 
to produce dimers from an atomic gas, already proposed 
in the case of bosons in 33| , with the disadvantage that 
the resulting molecules have a high center of mass kinetic 
energy, on the order of h? jmc? . These strong three-body 
losses prevent in present experiments to follow adiabati- 
cally the upper branch in Figure |21 from the left part of 
the Figure (a > 0) to the right part of the Figure (a < 0). 
In the unitary regime, our model predicts in the upper 
branch an energy per particle scaling as f? p^l^ jm larger 
than the ideal Fermi gas case, that is corresponding to 
an effective repulsion. A similar behavior was predicted 
for bosons by A. Leggett 0. 

Inspection of Figure |2] immediately inspires another 
way to produce a gas of ultracold molecules than the one 
relying on three-body inelastic collisions: one starts with 
a weakly attractive atomic Fermi gas, on the right of 
the Feshbach resonance, and one slowly crosses the res- 
onance from right to left in order to follow adiabatically 
the ground branch. Note that a similar mechanism has 
been developed in the case of an atomic Bose gas |26|. 
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FIG. 3: Pressure P of the gas in units of the Fermi pressure 
Pf = ti'kp jVa-K^m, as function of —l/kpa- Only the first 
two branches of the model are represented. 



Another important feature of the model is that noth- 
ing dramatic happens right on the Feshbach resonance, 
the mean energy per particle being simply proportional 
to the Fermi energy with some numerical factor (here 
3/40). As this mean energy is less than the ideal gas 
value l(2Jl, the gas experiences an effective attraction due 
to the atomic interactions. This universality and this ef- 
fective attraction appeared already in several approaches 
0, Jf ^J i ^ , and are confirmed by the experimental 
results 111- As the total energy remains positive, so 
will be the pressure and the compressibility of the gas at 
the Feshbach resonance: the Fermi gas will not experi- 
ence a collapse when kpa becomes large in absolute value 
and negative, in agreement with the recent experiments 
Q Q HI ^ IS and contrarily to what was feared a few 
years ago based on the calculation of the compressibility 
in the mean field approximation 27] . 

To check in a systematic way the stability of the gas 
in our model, we have calculated the pressure (see Fig- 
ure |2J) and the compressibility (not shown) of the gas for 
the lowest two branches for an arbitrary value of a. The 
compressibility and the pressure are positive everywhere. 
Moreover, the pressure on the ground branch is lower 
than the ideal Fermi gas pressure and tends exponen- 
tially to zero from above in the molecular limit a — 0^, 
revealing the absence of a non-zero scattering length of 
dimers in our model p^ . The decrease of the pressure 
across the Feshbach resonance from the a < part to the 
a > part is observable in a trap through a shrinking of 
the size of the cloud. 

In conclusion, the heuristic model that we have pre- 
sented contains all the essential features of a two spin- 
component Fermi gas for s-wave interactions with an 
arbitrary scattering length, such as the existence of a 
metastable atomic phase (upper branch) and a molecu- 
lar phase (lower branch) for a ^ O"*" and the continous 
connection between the molecular regime and the weakly 
attractive regime when a is varied from O"*" to 0~ across 
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the Feshbach resonance. 

Note: after the submission of this paper, both sce- 
narios for the production of a molecular condensate 
have been realized experimentally: the adiabatic scenario 
^9. .3Qj and the scenario assisted by three-body inelastic 
collisions [sills^l. 
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